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Abstract. 

Three dimensional (3D) galactic cosmic ray (GCR) anisotropy has been studied for 2006- 
2012. The GGR anisotropy, both in the ecliptic plane and in polar direction, were obtained 
based on the neutron monitors (NMs) and Nagoya muon telescopes (MT) data. We analyze 
two dimensional (2D) GCR anisotropy in the ecliptic plane and north-south anisotropy normal 
to the ecliptic plane. We reveal quasi-periodicities - the annual and 27-days waves in the GCR 
anisotropy in 2006-2012. We investigate the relationship of the 27-day variation of the GCR 
anisotropy in the ecliptic plane and in the polar direction with the parameters of solar activity 
and solar wind. 


1. Introduction 

Diurnal variation, being the measure of the galactic cosmic ray (GCR) anisotropy, arises from 
the Earth’s rotation as a ground-based detector’s viewing direction through the sky during each 
24-hours period. Three-dimensional (3D) anisotropy vector is determined by the distribution 
of the stream of cosmic rays in the 3D heliosphere. An average anisotropy vector has been 
explained based on the diffusion-convection theory of GCR modulation in the heliosphere [1]- 
[3]. That is a consequence of the equilibrium established between the radial convection by solar 
wind and the inward diffusion of GCR particles along the interplanetary magnetic field (IMF) 
owing to the radial gradient of GCR. 

Scientific staff of IZMIRAN’s cosmic ray laboratory 
(http;//hehos.izmiran.troitsk.ru/cosray/main.htm), have calculated the components Ar, Af, 
and A 0 of the 3D anisotropy by Global Spectrographic Method (GSM) Hj-[6] based 
on hourly data from all operating neutron monitors (NMs). Unfortunately, the 
derivation of the Aq component is possible with an accuracy up to constant, and so, 
a value of the north-south anisotropy Aq obtained by GSM method is not accurate 
(http://cr20.izmiran.rssi.ru/AnisotropyCR/index.php). That is, sorrowfully, a deficiency of the 
GSM related to the nature of NMs data. 

In general determination of the two dimensional (2D) anisotropy in the ecliptic plane is 
feasible based on establishing the radial Ar and tangential Af components by the harmonic 
analysis method for an individual detector (e.g., NM or Muon Telescope (MT)). Alania et ah, [7] 
demonstrated that the 2D GCR anisotropy calculated by the radial and tangential components 
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Figure 1. Daily data of the sunspot number SSN, radial Ar and tangential Af components 
and phase of the 2D GCR anisotropy for Oulu NM in 2006-2012. 
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Figure 2. Temporal changes of the (top) daily GG index for 2006-2012, (midle) annnal trend 
in GG index and (bottom) daily GG index with excluded annual trend. 


determined using GSM basically do not differ from the values of anisotropy found by harmonic 
analysis method for individual NM with cut-off rigidities < 5 GV near the solar minimum epoch. 

Swinson [2j first suggested that the north-south anisotropy could be related to the cosmic 
ray flow caused by a positive heliocentric radial density gradient Gr of cosmic rays and the By 
component of the interplanetary magnetic field (IMF) B, as expressed by the vector product, 
By X Gr. Since the average By and the Gr lie in the ecliptic plane, the direction of the vector 
product By x Gr (or the north-south anisotropy) is expected to be perpendicnlar to the ecliptic 
plane and to reverse direction with changes of the magnetic field direction. It is directed upward 
when By is positive (IMF away from the Snn-^ > 0) and vice versa when By is negative (IMF 
toward the Sun-A < 0). 

Papers m- m studied the cosmic ray anisotropy vector in three dimensions. They have 
determined north-sonth anisotropy nsing data from polar located NMs. They conclnded that 
a magnitude of the north-south anisotropy varies with ~ 11-years period and there is not any 
dependence on solar magnetic polarity. 

An alternative method to study north-south anisotropy is proposed by Mori and Nagashima 
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Figure 3. Daily GG index, By component and the sector structure of the IMF in 2008 


m- They have introduced an index GG calculated from Nagoya MT data as: 

GG = ^[(49iV-495) + (49iV-49F)] (1) 

The terms represent the counting rate differences for telescopes pointing in North, South, and 
East directions at 49*^ zenith angle. The GG index is free of noise in isotropic intensity caused 
by Forbush decreases, periodic variations, atmospheric temperature effects, and geomagnetic 
cut-offs. GG index mainly reflects change of difference between the intensity from north polar 
direction and that from north and parallel to the equatorial plane directions. Counting rate 
differences of the north-south and north-equatorial could not be precisely contained the same 
type of information. In spite, a GG is accepted by cosmic ray community as a good alternative 
index to study the north-south (or north-equatorial) asymmetry of the cosmic ray flux e.g., [14| . 

The present paper investigates the annual and 27-day waves of the GCR anisotropy based on 
ground observations of secondary cosmic rays by NMs and MTs near the recent solar minimum 
2006-2012. We investigate the relationship of the 27-day variation of the GCR anisotropy in the 
ecliptic plane and in the polar direction with parameters of solar activity and solar wind. 

2. Data 

Figure [U presents the temporal changes of the daily sunspot numbers (SSN), daily radial Ar 
and tangential Af components, and phase of the 2D GCR anisotropy obtained by Oulu NM 
for 2006-2012. Figure [1] shows that components of the daily ecliptic plane GCR anisotropy 
obtained from hourly data of NM experience large dispersions during the Sun’s rotation period. 
However, one can recognize that the average daily radial component Ar oscillate near zero, 
which is caused by a sign dependent drift effect in well established sector structure of the IMF; 
the tangential component Af (~ 0.4%) represents diffusive corotational anisotropy and phase 
(time of maximum intensity) indicates that GCR stream is directed averagely from ~ 18h. 

Figure [2] (top panel) presents daily data of the GG index for 2006-2012. Figure [5] shows 
that a clear annual trend is seen in changes of the GG index. 

Our aim is to study both annual wave and 27-day variation in GG index. We reveal annual 
trend by smoothing daily data over 27 days (Fig. [2] middle panel). Finally for studying the 
27-day variation of the GG index we exclude the annual trend and present detrended daily GG 
in the bottom panel of Fig. [2j 

As far the formulation of the north-south anisotropy is based on the drift model, it has 
been established that there is a good correlation between GG index and the polarity of the IMF 


































































Figure 4. Wavelet analysis of the daily solar wind velocity V (a), By component of the IMF 
(b), GCR intensity for Oulu NM (c) for 2007-2008. 



Figure 5. Wavelet analysis of the Ar (a) and Af (b) components of the 2D GCR anisotropy 
and daily GG index (c) for 2007-2008. Dashed line designates the period of 27 days. 


during each solar rotation m- Consequently the GG index being the measure of the north-south 
asymmetry is inversely related with By component of the IMF. It is also well pronounced in the 
recent solar minimum 23/24. As an example Fig. [3] presents daily GG index, By component 
and the sector structure of the IMF in 2008. 

3. 27-day wave in GCR anisotropy near the solar minimnm 23/24: 2007-2012 

The 27-day variation of the GCR anisotropy has been studied less intensively up to present, 
in general. It partially is connected with the complexity of the reliable revealing of the 27-day 
variation of the GCR anisotropy by means of small amplitudes of the diurnal variations of GCR 
(< 0.3%, measured by NMs), and with a large dispersion of the amplitudes of diurnal variation 
comparable with the accuracy of hourly data of NMs. In papers Alania et ah, 0 , m studied 




































































































































Figure 6. Wavelet analysis of the daily solar wind velocity V (a), By component of the IMF 
(b), GCR intensity for Oulu NM (c) for 2009-2010. 



Figure 7. Wavelet analysis of the Ar (a) and Af (b) components of the 2D GCR anisotropy 
and daily GG index (c) for 2009-2010. Dashed line designates the period of 27 days. 


the 27-day variation of the 2D GCR anisotropy in the ecliptic plane. They demonstrated that 
the average amplitude in the minimum epoch of solar activity is polarity dependent, as it is 
expected from the drift theory. 

The 27-day variation of the north-south anisotropy was studied in the series of papers 
by Swinson and coauthors |20j-|22j. They showed that 27-day variation of the north-south 
anisotropy is correlated with solar activity and this correlation is not clearly dependent upon 
solar magnetic polarity. 

The recent solar minimum provided an unique opportunity to study recurrent variations under 































































































































Figure 8. Wavelet analysis of the daily solar wind velocity V (a), By component of the IMF 
(b), GCR intensity for Oulu NM (c) for 2011-2012. 
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Figure 9. Wavelet analysis of the Ar (a) and Af (b) components of the 2D GCR anisotropy 
and daily GG index (c) for 2011-2012. Dashed line designates the period of 27 days. 


relatively stable low solar activity conditions. Recurrent variations connected with corotating 
structures (~ 27 days) are clearly established in all solar wind and interplanetary parameters. 
Consequently the 27-day recurrent variations of cosmic ray intensity were clearly seen in a variety 
of cosmic ray counts detected by neutron monitors (e.g., m, [E]) and space probes (e.g., [IZj). 
Recently, Yeeram et ah, |23] studied recurrent ’trains’ (trend) of enhanced GCR anisotropy 
under influence of corotating solar wind structures near the recent solar minimum. 

To reveal the quasi-periodic variation of the 2D GCR anisotropy in ecliptic plane 
and north-south asymmetry of GCR stream by means of GG index near the recent 

















































































































solar minimum 2007-2012, we employ the wavelet time-frequency spectrum technique 
developed by Torrence and Compo [23]. The wavelet software is available at the website 
Ihttp://paos.colorado.edu/research/wavelets/software.htmll. In our calculation we used the 
Morlet wavelet mother function. 

We performed wavelet analysis for the daily solar wind velocity V, By component of the 
IMF, GCR intensity for Oulu NM, GG index, Ar and Af components of the 2D OCR anisotropy 
during 2007-2012 considering as a sampling time interval equaling two years (2007-2012). In 
our case a time interval of two years gives a good enough statistics as far our aim is to reveal 
recurrences < 27-day period, which is ~ 4% from whole sampling period of 26 solar rotations 
(two years). 

Results of calculation for V, By and GCR intensity are presented in Fig. Hjfor 2007-2008, 
Fig. El for 2009-2010, Fig. [5| for 2011-2012, respectively. Wavelet analysis of Ar and Af 
components of the 2D GCR anisotropy for Oulu NM and GG index are presented in Figs. [5|for 
2007-2008, Fig. [7|for 2009-2010, Fig. [9] for 2011-2012, respectively. Figures [5l [3 [9] panel (c) 
present very clear quasi-periodic changes in GG index related to the Sun’s rotation (~ 27 days) 
for almost whole time interval 2007-2012. Similar quasi periodic character is clearly visible in V 
and By component of the IMF. 27-day variation of GCR intensity is well established at the end 
of 2007 and in 2008. Although recurrent variations connected with corotating structures (~ 27 
days) are clearly established in almost all solar wind parameters, 2D GCR anisotropy shows a 
weak 27-day variation, only in some periods. This is connected with large dispersion of daily 
Ar and Af components of the 2D GCR anisotropy. 

4. Annual waves in GG index near the solar minimum 23/24: 2006-2012 

We have also studied a feature of the north-south asymmetry based on the GG index depending 
on the position of the Earth (±7.5^ in September and March) during its annual motion around 
the Sun. For this purpose we calculated 3 month average of the GG index for each year during 
2006-2012. In each year 1st point is the average of November-December-January, 2nd point 
is the average of February-March-April, 3rd point is the average of May-June-July and 4th 
point is the average of August-September-October. Results of calculations for superimposed 
GG index are presented in Fig. [TUJ Figure [TUI shows that GG index exhibits clear annual wave 
with evident maximum phase near the February-March-April period. We suppose that it is an 
indication of the asymmetry of a heliolatitudinal distribution of solar activity having minimum 
in south hemisphere near 7.5®; being some source of north-south asymmetry of the GCR flux 
in 2007-2012. Unfortunately, we could not make a precise quantitative estimation due to some 
uncertainty character of GG data used in this paper. At the same time, one can suppose that 
this finding is an indication of an existence of north-south asymmetry of the GCR flux. 

5. Conclusions 

(i) We study recurrent variations of the GCR anisotropy in the ecliptic plane and in the polar 
direction connected with corotating structures observed in the heliosphere near the recent 
solar minimum in the period 2007-2012. 

(ii) 2D GCR anisotropy generally does not show a clear evidence of the 27-day variation, but 
in some periods weak recurrent character is visible. This is partially connected with large 
dispersion of daily Ar and Af components of the 2D GCR anisotropy. 

(iii) Using wavelet time-frequency method we reveal clear 27-day waves in the GG index, being 
the measure of the north-south anisotropy, for almost whole analyzed period 2007-2012. 
Similar recurrent properties are observed in solar wind velocity and interplanetary magnetic 
field. 
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Figure 10. Superposition of the 3 month average GG index for each year during 2006-2012 


(iv) We show that GG index exhibits a clear annual wave with an evident phase of maximum 
near the February-March-April period. One can suppose that possibly it is an indication of 
the asymmetry of a heliolatitudinal distribution of solar activity having minimum in south 
hemisphere near ~ 7.5*^; this heliolatitudinal asymmetry can be considered as a some source 
of north-south asymmetry of the OCR flux in 2007-2012. 
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